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[57] ABSTRACT 

A radiographic scanner (A) generates a high energy 
image representation which is stored in a high energy 
image matrix (V) and a low energy image representa- 
tion which is stored in a low energy image memory (U). 
A pair of filter functions selecting circuits (C) select a 
first or soft tissue specific filter function and second or 
bone specific filter function, respectively. The soft tis- 
sue filter function selecting circuit selects and adjusts 
the soft tissue filter function in accordance with the 
pixel value of the low energy image representation for 
each corresponding pair of pixel values. Convolvers 
(44, 46) convolve pixel values from the high and low 
energy image representations with the selected and 
adjusted filter functions. A soft tissue transform func- 
tion (48) transforms the filtered high and low energy 
image representations into a soft tissue or other material 
specific image representation (42). The other filter se- 
lecting and adjusting circuit selects and adjusts the bone 
specific filter functions which are convolved with the 
high and low energy image representations by convolv- 
ers (54, 56). A bone specific transform function (58) 
transforms the filtered high and low energy image rep- 
resentations into a bone basis image. 

21 Claims, 3 Drawing Sheets 
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taken simultaneously with back to back detectors, one 

ADAPTIVE FILTER FOR DUAL ENERGY of which was sensitive primarily to the lower energy 

RADIOGRAPHIC IMAGING x-rays and the other was sensitive primarily to the 

higher energy x-rays. 

TECHNICAL FIELD 5 In the prior art dual energy digital radiography, both 
This invention relates to the art of diagnostic imag- *"gh energy and low energy electronic image represen- 
ing. It finds particular application in conjunction with tations were derived. Each representation commonly 
dual energy, medical diagnostic digital radiography and included a rectangular array of pixels, each pixel having 
will be described with particular reference thereto. a pixel value indicative of the degree of radiation atten- 
However, it is to be appreciated that the invention may 10 uation or transmissivity through a corresponding path 
find further application in conjunction with other imag- of the region of interest. The pixel value variations were 
ing techniques such as computed tomography and the commonly represented by gray scale variations in man- 
like, readable images. Utilizing the transform functions set 
DAP^roniixm AT)T forth in "Generalized Image Combinations in Dual 
BACKGROUND AKI is Kvp Digita i Radiography," L. A. Lehmann, et aL, 
In a conventional film radiography system, an x-ray Medical Physics, Vol 8, No. 5, pages 659-667, Septem- 
source directed a divergent beam of x-rays through a ber/October 1981, the high and low energy specific 
patient After passing through the patient, the x-ray images were transformed into one or more material 
beam irradiated a phosphorescent screen as well as a specific basis image representations. The Lehmann, et 
light and x-ray sensitive film positioned adjacent the 20 al. transform operated on the pixels of the high and low 
screen. As the radiation passed through the patient, it energy images that correspond to the same volumetric 
was attenuated in accordance with the tissue through subregion of the region of interest to produce a corre- 
which it passed to produce a shadow image on the film. sponding pixel value of the material specific image. 
The radiologist examined the gray scale, i.e. the light Mo8t commonly, two basis images were generated — 
and dark regions, of the film directly for diagnostic « Qne for of soft ^ ^ othef for bonc of 

purposes. calcium. For calibration purposes,- plexiglass was uti- 

In digital radiography, the gray scale of each pixel or lized to approximatc water ^ soft tissuc md 

incremental region of the Hnagc : was represented by an Wfl3 utilkcd approximate bone ^ 0tner 

electronic, digital value. These Oimtal values were pro- ™ ch ^ ^ by ^ subtraction of 

cessed by various data processing and image enhance- 30 * 

ment techniques to improve the diagnostic value of the ^j^d cu «s* un* S ca, uavc wcu 

dl c 0ne . of the P r ° blcms ™* the Lehmann et al. trans- 

ri£ the electronic digital values f onn f ° K n wa ? that lt ^f^f d n0lse , ° f the image* 

In one prior art digital radiographic technique, the 35 |f ^^ fc nn ?» e ? lc ! wer "W^-to-noisc ratio 

digital values were derived directly from photographic clther . or low energy images. The 

film using a photodensiometer, video camera, or the ™* om noise degradation was particularly bad m the 

like. In another prior art technique, a phosphor with a *™y attenuation or thick portions of the spine 

long term memory temporarily recorded the shadow- wluch were represented by relatively low pixel values, 

graphic image. A photo-sensitive apparatus converted 40 Various techniques have been developed for operat- 

the phosphor luminescence into corresponding digital m 8 on the soft tissue and bone specific basis images to 

values on a pixel by pixel basis. reduce the effect of random noise. One technique was to 

The intensity of x-ray beams traversing the patient convolve each pixel of the basis image with a smoothing 

have also been converted directly into electronic sig- or fflt « function that averaged each pixel value with a 

nals. An array of electronic radiation detectors was 45 percentage of the average value of the surrounding 

disposed opposite the patient to receive the x-rays pass- P«els. One drawback of this technique was that it re- 

ing therethrough. In one technique, often denoted as moved the apparent random noise in the dense bone 

"scan or slit projection radiography", a thin fan beam of regions of the image by blurring the entire image, 

radiation passed through a narrow plane of the patient To avoid blurring all regions of the basis image, oth- 

and impacted a linear array of x-ray detectors. The 50 ers have provided a region specific filter function. Ran- 

radiation fan and the detector array were moved trans- dom noise commonly appears as a speck of a different 

verse to the fan beam plane to scan a selected rectangu- intensity from the surrounding areas, i.e. a pixel value 

lar region of the patient. that is much larger or smaller than the values of adja- 

In another direct detection technique, an "area" beam cent pixels. Specifically, the deviation between a pixel 

of radiation was directed through the entire rectangular 55 value and each of its surrounding neighbors was dcter- 

region of interest simultaneously to impact upon a large mined. If greater than a preselected deviation was ob- 

rectangular detector array. The intensities detected by served, the blurring filter function was applied to that 

the detectors were digitized for appropriate data pro- pixel of the basis image. If less than the preselected 

cessing and display on a video monitor or the like. deviation was observed, no filtering was applied to that 

The material or body tissue can be characterized by 60 pixel, 
the difference in its attenuation of a high energy and a One disadvantage of this selective filtering technique 
low energy x-ray beam. To improve the diagnostic is that the deviation between neighboring pixel values 
value of radiographic images, images of the region of was not always an indication of noise. Rather, the devia- 
interest have been constructed from both high and low tion may have been the result of a tissue or bone bound- 
energy radiation. In KV switching, dual energy radiog- 65 ary, a thin blood vessel, or the like. This led to an aver- 
raphy, two time sequential images were taken, one with aging which tended to blur edges and obscure fine de- 
higher energy x-rays and the other with lower energy tails of the filtered basis image. Other filtering and data 
x-rays. In dual detector scanners, both images were processing techniques, such as edge enhancement tech- 
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niques, other types of smoothing or filtering, and the vided. First and second image memories are provided 

like have also been performed on the resultant basis for storing first and second image representations of the 

images. same region of interest of the patient. A filter means 

In accordance with the present invention, there is separately filters pixel values of the stored first and 

provided a new and improved image enhancement 5 second image representations with a first filter function 

method and apparatus which overcomes the above to generate first and second filtered image representa- 

referenced problems and others. tions, respectively. A transform means transforms the 

SUMMARY OF THE INVENTION t^^t^ ** "? T"* T 

sentations into a composite image. A composite image 

In accordance with one aspect of the present inven- 10 memory is connected with the transform means for 

tion, a method of diagnostic imaging is provided. First storing the composite image. 

and second image representations of the same region of In accordance with another aspect of the present 
a subject are generated. In the preferred embodiment, invention, an apparatus for diagnostic imaging is pro- 
one of the images is a higher energy x-ray image repre- vided. An image memory stores pixel values for each of 
sentation and the other is a lower energy x-ray image 15 a plurality of pixels. A filtering means operates on the 
representation. *Prior to operating on the first and sec- pixel values of each selected pixel and pixels thereadja- 
ond image representations with a transform function, cent with a selectively alterable filter function. A filter 
such as the Lehmann, et al. transform function dis- function altering means alters characteristics of the 
cussed above, the first and second image representations filter function in accordance with the pixel value of 
are each filtered. The first and second filtered image 20 each pixel to be filtered. In this manner, the filter func- 
representations are then transformed with one or more tion with which each pixel is filtered varies with its 
transform functions to produce a corresponding number pixel value. 

of composite images. If the Lehmann, et al. transform A primary advantage of the present invention is that 

functions are selected, two composite images com- it enhances electronic images and reduces noise degra- 

monly known as a soft tissue basis image and bone tissue 25 dation. 

basis image are produced. Other transformations of the Another advantage of the present invention is that it 

filtered first and second image representations, such as adaptively filters the noisier regions of an image more 

simple subtraction, weighted averaging, approximations than more noise free regions. 

of the Lehmann, et al. equations, and the like may also Noisier thick bone representing regions are filtered to 

be utilized. 30 remove noise while substantially noise-free regions rep- 

In accordance with another aspect of the present resenting the lungs are not blurred or filtered to assure 

invention, a method of diagnostic imaging is provided. their diagnostic integrity. 

A generated image representation is defined by a plural- Still further advantages of the present invention will 
ity of pixels, each pixel having a pixel value. For one of become apparent to others upon reading and under- 
the pixels, its pixel value is compared with a preselected 35 standing the following detailed description of the pre- 
criteria, such as a preselected pixel value amplitude. ferred embodiments. 
One of a plurality of preselected filter functions is se- 
lected in accordance with the comparison of the pixel BRJEF DESCRIPTION OF THE DRAWINGS 
value and the preselected criteria. The pixel values of The invention may take form in various steps and 
said one pixel and adjacent pixels are filtered with the 40 arrangements of steps or in various components and 
selected filter function to create a filtered pixel value for arrangements of components. The drawings are only 
a corresponding pixel of a filtered image representation. for purposes of illustrating a preferred embodiment and 
The comparing, filter selecting, and filtering steps are are not to be construed as limiting the invention, 
repeated for each pixel of the image representation to FIG. 1 is a diagrammatic view of a system for gener- 
produce the corresponding pixel values for the filtered 45 ating soft tissue and bone specific images from dual 
image representation. In this manner, the filter function energy radiographic data in accordance with the pres- 
for each pixel is selected in accordance with that pixel's ent invention; 

pixel value. This method finds particular application in FIGS. 2A and 2B illustrate a preferred filter function 

conjunction with digital x-ray radiography techniques alteration scheme for the soft tissue image of FIG. 1; 

in which the signal-to-noise ratio decreases as the mag- 50 FIGS. 3A and 3B illustrate a preferred bone tissue 

nitude of the pixel value decreases. filter function altering scheme; and, 

In accordance with another aspect of the present FIG. 4 is an alternate embodiment to the system of 

invention, a method of medicaT diagnostic imaging is FIG. 1 for producing a combined soft tissue and bone or 

provided. First and second image representations of a "normal" image, 
common region of a patient are generated. Both the first 55 

and second image representations are operated on with DETAILED DESCRIPTION OF THE 
(1) at first transform function and (2) a first filter func- PREFERRED EMBODIMENTS 
tion to create a filtered first composite image, preferably With reference to FIG. 1, a diagnostic scanner A 
a soft tissue basis image. Both the first and second image such as a digital, x-ray slit scanner, a CT scanner, a 
representations are also operated upon with (1) a second 60 magnetic resonance imager, or the like generates a pair 
transform function and (2) a second filtering function of image representations of the same region of interest 
which is different from the first filtering function to of a subject, e.g. a low energy image U and a high en- 
create a filtered second composite image, preferably a ergy image V. An adjustable filter and transform means 
bone basis image. In-this manner, different filter func- B converts the high and low energy images into one or 
tiona.are.utilized in the generation of the bone and soft 65 more material specific or other composite images. A 
tissue indicative images. filter altering means C monitors at least one of the high 
In accordance with yet another aspect of the present and low energy images to determine the relative magni- 
invention, an apparatus for diagnostic imaging is pro- tude of the signal-to-noise ratio at each image pixel and 
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alters the function applied in the filter and transform The high and low energy image representations are 

means B accordingly. The lower the signal-to-noise each defined by a rectangular array of pixels, corre- 

ratio is in a subregion, the more the pixel values repre- sponding to a rectangular array of subregions of the 

senting the subregion are filtered or smoothed scanned region of interest. The digital pixel values cor- 

The scanner A illustrated in FIG. 1 for purposes of 5 responding to each pixel of the images or subregions of 

explanation is a slit scan projection type digital radiog- the region of interest are stored in corresponding mem- 

raphy system. It includes an x-ray source 10, such as an orv elements of low and high energy image memories U 

x-ray tube, which generates a polyenergetic fan beam 12 and V. In the illustrated embodiment, one pixel in each 

of x-rays having at least two energy levels. An entrance column of the images corresponds to a corresponding 

collimator 14 and an exit collimator 16 collimate the 10 detector in the detector array. Each time the detector 

x-ray beam passing through the patient to be predefined ^ swept m incremental transverse distance, the 

fan shaped plane. A first or low energy detector array data for ^ next adjacent pixel commn is generated. 

18 includes a first linear array of detectors for detecting ^ detector amys m back t0 back> ea ch pixel 

the low energy x-rays of the x-ray beam A second or . ^ low ^ reprcscntation ^ a corrc . 

high energy x-ray detector array 20 inc hides a hnear 15 s ding ^ ^ ^ ^ energy ^ representation 

array of x-ray detectors disposed d^ectly behmd the ^ ^ ^ ^ subregion of th ; regi on of 

low energy detector array to detect the high energy ^ ^ preferred embodiment, subregions of the 

x-ray photons. subject which attenuate radiation the most, such as 

The stopping power of the low enei -gy detectors^ ^ c ^ ^ scintillation tals t0 te 

such that pr^ardy low energy pho^ 20 hence J ^ 

transformed to electnca^ signak, whUe most .fogh en- pixel P valucs . other numbering schemes for 

ergy photons pass *^ ou f\^ the pixel values may, of course, be utilized. The values 
tectors are configured of a denser higher atomic number r . : , ■ « -*t. *u 
material with more stopping power to stop and convert ™3\ var * l°gnthmically or with other prese- 
the high energy photons into electrical signals indica- 25 ^ ^f^T* ™ y T™*' * . t A 
tive of the mtensity of high energy radiation impinging ™ e adjustable filter and transform means^B mcludes 
thereon. One detector for each of the high and low * firs f m ^ W fdt f transform means 
energy arrays are disposed in alignment back to back 40 whlch fflter _ s transfonns the high and low energy 
with the x-ray source such that a corresponding pair of ^ mt0 a finit or tlssuc s t tored m 
the detectors detect high and low energy radiation 30 memory 42. A fu^t filter means convolves the low and 
which has traversed the same line or subregion of the ^ enQT ^ representations separately pixel by 
patient. An x-ray absorptive material 22, such as a sheet with an adaptive first filter function that varies 
of copper, absorbs middle energy x-ray photons that from pixel to pixel to produce filtered low and high 
have sufficient energy to pass through the low energy energy image representations. A first convolver 44 op- 
detectors. The thickness of the copper sheet is selected 35 crates on the low energy image pixel values with the 
to absorb photons whose energy is below a preselected filter function to produce the filtered low energy 
high energy level. In this manner, relatively distinct and image representation. A second convolver 46 operates 
separated high and low energy photon levels are de- oa eac ^ of the high energy image pixel values with the 
tected. same first filter function as the corresponding low en- 

A mechanical drive 24 causes the high and low en- 40 ergy image pixel to generate the filtered high energy 
ergy detector arrays and the entrance and exit collima- image representation. That is, the low energy image 
tors to be swept transverse to the plane of the beam such pixel value for the ith column and the jth row, U(i j), 
that a rectangular region of interest of the patient is and the high energy image pixel value for the same 
irradiated. An encoder 26 changes the address for the column and row, V(iJ), are each convolved with the 
low and high energy images U and V each time the 45 same filter function K(iJ) to generate filtered values 
detector array moves an incremental distance. In this U'(ij) and V'(ij) corresponding to the (i j)th pixel. Op- 
manner, sweeping the detector array provides analo- tionally, a filtered image memory means may be pro- 
gous data to a pair of rectangular high and low energy vided for storing the filtered high and low energy im- 
detector grids. Optionally, such detector grids may be ages from the convolvers 44 and 46. 
provided with the elimination of the mechanical drive 50 A first transform means 48 operates on the filtered 
24 and encoder 26 such that all pixels of the high and low and high energy image representations U' and V 
low energy images are collected simultaneously rather with the first transform F to generate a selected material 
than column by column. specific image I. That is, the material specific image I is 

Each detector of the high and low energy detector a function F of the filtered high and low energy images, 

arrays produces an analog electronic output signal 55 Tj' md y : 
which varies in proportion with the intensity of x-rays 

which it receives. In the preferred embodiment, each l=F{U.V) (i) 
detector includes a scintillation crystal which scintil- 
lates or emits light in proportion to the intensity of the Various functions F may be selected. As one exam- 
received radiation. A photodiode produces an analog 60 pie, each pixel value of the material specific image I . 
output signal whose amplitude varies in accordance may be the difference of the corresponding pixels of the 
with the intensity of the light from the scintillation high and low energy images: 
crystal. Detectors which convert the radiation directly 

into electrical signals, CCD devices, and other known AW) = (2) 
radiation detectors are also contemplated. Analog to 65 

digital converters 28 convert the analog intensity sig- As yet another example, the material specific image 

nals into digital representations of the radiation inten- may be other linear combinations of the filtered high 

sity. and low energy images, such as a weighted averaging: 
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value increases, hence, the radiation attenuation de- 
ran _ m van + n van (3) creases, the signal to noise level ratio improves and the 

"> + n * amount of filtering is reduced. 

In the preferred embodiment, the Miter function K is 
where m and n are weighting factors. 5 defined by a 3 X3 kernal matrix: 

In the preferred embodiment, a transform analogous 
to the transform set forth by Lehmann, et al. is utilized. 
Each image of the material selected image is related 
logrithmically to the corresponding pixels of the fil- 
tered high and low energy images: io 



B B 
3 A 
B B 



(5) 



m - *iiog ir«/> + Ajlog vu + < 4 > convolving a given pixel value U(ij) with this filter 

kilos 1 V(U> + *4og 2 r(v) + function, replaces the given pixel value with the sum of 

Aaiog cr(v) log V{ij), I 5 A times ^ S* vcn P^ 1 value P lus B of tne 

eight immediately adjacent pixel values, i.e. the filtered 

where ki, k 2 , k 3 , k* and k 5 are constants which are U ' is defined by: 
determined for the given scanner A and for a selected 

material to which the image is to be selective in an m AU ^ + ^ _ ij _ i) + BU(i - \$ + (6) 

initial calibration proems. Due to the hardware con- ™ BU{i _ , + 1} + . 

stramts, the k constants will be different for each se- nim ^ . . 1W , , ' 

lected scanning apparatus. Moreover, different con- - t ! V7 J> a + + t l <» + + + !) 

stants will be selecSd for a soft tissue selective image, a JA"S B ?* ^T* ? W1 * 

bone selective image, and other material selective im- ^ "mptude of the given pixel value. The preferred 

ages. Plexiglass may be utilized for an initial calibration 25 embodiments o £ and B for soft tissue are defined by 

to select the k constants for a soft tissue specific image " ™ of I* : ^ c , 

and aluminum may be utilized as a calibration phantom u Wl . th f^ lfic . reference to FIGS. 2A and 2B, when 

to select the calibration constants k for a bone specific ™ pwel valueis below a preselected minimum thresh- 

image. old 70, about 2% of the pixel value dynamic range in the 

A second material specific filter and transform means 30 preferred embodiment, the region is assumed to be very 
50 filters and transforms the high and low energy im- nois y* Accordingly, A and B are each selected to be 
ages into a second, preferably a bone specific image x/ ?< ***** to » ^ P** 61 value is replaced with the average 
stored in memory 52. A pair of convolvers 54 and 56 of itself and its eight surrounding pixel values. When the 
operate pixel by pixel on corresponding pixels of the P** el value is above a preselected high threshold 72, 
high and low energy image representations with a sec- 35 preferably about 20% of the pixel value dynamic range, 
ond adaptive filter function K2. The convolvers 54 and toe P** 6 * vame » unchanged. That is, A = 1 and B =0. 
56 perform a similar convolution as convolvers 44 and ^ n between the high and low threshold values, the 
46 except with a different filter function that is uniquely relative magnitudes of A and B vary linearly with the 
selected for the material specific image to be generated. P**^ levels. In the embodiment illustrated in FIGS. 2A 
A second transform means 58 operates on the filtered 40 and 2B, the value of A varies linearly between 1/9 and 
high and low energy images with a second transform G one; the value of B varies linearly between 1/9 and zero, 
to transform them into a second material selective im- With reference to FIGS. 3A and 3B, a different filter 
age, preferably a bone selective image I. A display function is preferably selected for the soft tissue and 
means 60 selectively displays the soft tissue and bone bone specific images. Rather than the relationship illus- 
specific images from memories 42 and 52. As discussed 45 trated in FIGS. 2A and 2B for soft tissue, the preferred 
above, the function G may be a difference, weighted values of A and B for a bone specific image are illus- 
average, logrithmic combination, or other linear or trated in FIGS. 3A and 3B. For low pixel values or high 
non-linear combination of the corresponding pixels of x-ray attenuation below a low threshold 74, each pixel is 
the filtered high and low energy image representations. again averaged evenly with its surrounding neighbors. 
Preferably, the second transform G implements Equa- 50 Above a preselected high threshold level 76, each pixel 
tion (4) except that the values of the constants k are continues to be averaged with its surrounding neigh- 
selected specifically for bone or aluminum rather than bors. However, the value of A is selected at i and the 
soft tissue or plexiglass, value of B is selected at 1/16. Between the high and low 

Where appropriate, additional material specific fil- threshold values, the value of A varies between 1/9 and 

ter/transform means may be provided for filtering and 55 J linearly and the value of B varies linearly between 1/9 

transforming the high and low energy images into im- and 1/16. 

ages which are specific to other materials. Although the Other relationships between A and B may be se- 

filter convolution and the transform are shown as sepa- lected. However, it is preferred that: 
rate operations, it is to be appreciated that some filter 

and transform functions may be combined into a single 60 v4+8£«coastant (7) 
adaptive function to expedite data processing. 

The filter altering means C projects the signal-to- In the illustrated embodiment, the constant is one. 

noise ratio at each pixel of the image in accordance with Other relationships between A and B may clearly be 

the level of the pixel value and alters the filter function selected, for example, their values may change expo- 

or the transform accordingly. The lower the pixel 65 nentiaUy, logrithmically, or along other complex func- 

value, i.e. the greater the radiation attenuation, the tions. For example, A and B may vary with a generally 

greater the potential that the pixel value will include S-shaped curve in which A asymptotically approaches 

significant noise or be in a noisy region. As the pixel 1/9 at one extreme and 1 for soft tissue and 1 for bone 



05/11/2004, EAST Version: 1.4.1 



4,792,900 



at the other extreme; and B asymptotically approaches 
1/9 at one extreme and zero for soft tissue and 1/16 of 
bone at the other extreme. This function may vary con- 
tinuously over the full dynamic range of pixel levels 
without cut off thresholds. 

Other filter functions or kernals K may also be se- 
lected For example, the kernal may be other sizes, such 
as a 5X5 matrix which averages each pixel value with 
some fraction of the eight immediately surrounding 
pixel values and the sixteen next most immediately sur- 
rounding pixel values. Moreover, as the kernal matrix 
becomes larger, the pixel values within each surround- 
ing ring may vary generally in accordance with its 
radial distance from the center rather than solely as a 
function of the square ring. Other convolution and filter 
functions may optionally be selected. 

With reference again to FIG. 1, the preferred noise 
level detennining and filter altering means C includes 
separate circuitry for each of the material specific im- 



10 



10 



15 



tors the amplitude of pixel values of one of the images, 
preferably one of the filtered tow energy images. Vari- 
ous transform functions may be selected. The image 
created by the transform means 130 is stored in a com- 
posite image memory 134. For example, the transform 
means 130 may sum corresponding pixel values of all 
four filtered images when the corresponding pixel value 
of the monitored image is in a preselected amplitude 
range, such as a central range. For low amplitude 
ranges, the corresponding pixels of the two filtered high 
energy images may be added. In between the central 
and low ranges, a weighted averaging may be per- 
formed. Analogously, in a high pixel amplitude range, 
the corresponding pixels of the low energy image might 
be added. Between the middle and high pixel value 
ranges, a ramp function or the like might be selected as 
a basis for performing a weighted averaging. Other 
combination equations are, of course, contemplated. 
The invention has been described with reference to 
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filter function of FIGS. 2 and 3, to be selected for each 
material specific image. However, because each of the 
filter altering means C is of the same construction, it is 
to be appreciated that the following discussion applies 
equally to both. 25 

Where appropriate, an address scaler means 80 may 
adjust the low energy pixel values U(iJ) to correspond 
to the addresses of a kernal or filter function selection 
means 82. Preferably, the low energy image representa- 
tion is used to select the filter function for the corre- 30 
spending pixels of both the low and high energy images. 
In a preferred embodiment, the filter function selection 
means includes a first look-up table 84 which is ad- 
dressed by the first pixel value U(ij) to look up the 
corresponding A value and a second look-up table 86 35 
which is addressed by the pixel value to look up the 
corresponding B value. A kernal matrix means 88 as- 
sembles the kernal matrix or filter function K from the 
values retrieved from the look-up tables 84 and 86. If 
larger kernal matrixes are used, additional look-up ta- 40 
bles may be provided for selecting the additional matrix 
values. Optionally, the filter function selection means 82 
may be an arithmetic unit which calculates the A and B 
values mathematically from each pixel value U(i j). 

With reference to FIG. 4* the first and second or high 45 
and low energy images may be recorded at different 
times. An image digitizer 100 may convert a film, long 
term memory phosphor, or image intensifier output 
image into pixel values for corresponding pixels of a 
first or low energy image for storage in a first image 50 
memory means 102. Subsequently, a second film, phos- 
phor, or image intensifier image taken through the same 
region of the patient may be digitized by the digitizer 
for storage in a second image memory 104. If the first 
and second film or phosphor images were taken with 55 
different energies of radiation, the first and second digi- 
tal images will represent high and low energy images of 
the same region. Appropriate indexing to align the im- 
ages may be required. 

In accordance with the level of pixel values in the 60 
low energy or first image, filter function selecting 
means C select the appropriate filter functions for con- 
volvers 110, 112, 114, and 116 which generate filtered 
images for storage in filtered image memories 120, 122, 
124, and 126. 65 

An adaptive transform means 130 combines the fil- 
tered images in accordance with the level of pixel val- 
ues therein. A pixel value monitoring means 132 moni- 



and alterations will occur to others upon reading and 
understanding the preceding specification. It is intended 
that the invention be construed as including all such 
modifications and alterations insofar as they come 
within the scope of the appended claims or the equiva- 
lents thereof. 

Having thus described the preferred embodiments, 
the invention is now claimed to be; 

1. A medical diagnostic imaging apparatus for pro- 
ducing soft tissue and bone specific basis images from 
dual energy radiation data, the apparatus comprising: 
a radiation source for producing at least higher en- 
ergy and lower energy radiation; 
a radiation detection means for separately detecting 
the higher and lower energy radiation traversing a 
selected region of interest and for producing digital 
electronic values indicative of the intensity of 
higher and lower energy radiation received 
thereby; 

a higher energy image memory means operatively 
connected with the detection means for storing the 
higher energy radiation intensity indicative digital 
electronic pixel values for each pixel for a higher 
energy image representation of the region of inter- 
est; 

a lower energy image memory means operatively 
connected with the detection means for storing the 
lower energy radiation intensity indicative digital 
electronic pixel values for each pixel for a lower 
energy image representation of the region of inter- 
est; 

a first convolver means for convolving a first filter 
function pixel by pixel with the digital values of the 
higher energy image representation and for con- 
volving the first filter function pixel by pixel with 
the digital values of the lower energy image repre- 
sentation; 

a first filter function altering means for correspond- 
ingly adjusting a characteristic of the first filter 
function in accordance with the digital value of 
each pixel of at least one of the higher and lower 
energy image representations; 

a first transform means for transforming convolved 
values of corresponding pixels of the higher and 
lower energy representations into a digital gray 
scale value of a corresponding pixel of a soft tissue 
specific basis image; 



05/11/2004, EAST Version: 1.4.1 



a soft tissue image memory means operatively con- 
nected with the first transform means for storing 
the gray scale values of the soft tissue basis image; 

a second convolver means for convolving a second 
filter function with the digital values of the higher 5 
energy image representation and convolving the 
second filter function with the digital values of the 
low energy image representation; 

a second filter function altering means for corre- 
spondingly adjusting a characteristic of the second 10 
filter function in accordance with the digital value 
of each pixel of at least one of the high and low 
energy image representations; and, 

a second transform means for transforming con- 
volved values of corresponding pixels of the higher 15 
and lower energy image representations from the 
second convolver means into a digital gray scale 
value of a corresponding pixel of a bone specific 
basis image. 

2. An apparatus for producing diagnostic images, the 20 
apparatus comprising: 

a first image memory means for storing an electronic 
pixel value for each pixel of a first image represen- 
tation; 

a second image memory means for storing an elec- 
tronic pixel value for each pixel of a second image 
representation; 

a first filter means for operating on each pixel value of 
the first image representation from the first mem- 3Q 
ory means with a first filter function and for operat- 
ing on each pixel value of the second image repre* 
sentation from the second memory means with the 
first filter function; 

a first transform means for transforming filtered pixel 35 
values from the first and second memory means 
into a first composite image; and, 

a first composite image memory means for storing the 
first composite specific image. 

3. The apparatus as set forth in claim 2 further includ- ^ 
ing a first filter function altering means for selectively 
altering the first filter function in accordance with a 
statistical probability that the pixel to be filtered is free 
from noise degradation. 

4. The apparatus as set forth in claim 3 wherein the 45 
filter function altering means includes at least one look- 
up table which is operatively connected with at least 
one of the first and second image memories to be ad- 
dressed by the pixel value of at least one of each corre- 
sponding pair of pixels to be filtered, the look-up table 50 
retrieving kernal matrix values of a filter function in 
accordance with the addressing pixel value and wherein 
the first filter means includes a convolving means for 
convolving the first and second image memory pixel 
values with the kernal matrix. 55 

5. The apparatus as set forth in claim 2 further com- 
prising: 

a second filter means for operating on each pixel 
value of the first image representation from the first 
memory means with a second filter function and 60 
for operation on each pixel value of the second 
image representation from the second memory 
means with the second filter function; 

a second transform means for transforming filtered 
pixel values from the first and second memory 65 
means into a second composite image; and, 

a second composite image memory means for storing 
the second composite image. 
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6. The apparatus as set forth in claim 5 wherein the 
first and second image representations are radiation 
energy specific and the first and second composite im- 
ages are irradiated material specific. 

7. The apparatus as set forth in claim 6 further includ- 
ing a filter function altering means for selectively alter- 
ing the first and second filter functions, the filter func- 
tion altering means being operatively connected with 
one of the first and second image memories for altering 
the first or second filter function in accordance with 
each pixel value to be filtered by the filtering means. 

8. An apparatus for diagnostic imaging comprising: 
a first image memory which stores a plurality of pixel 

values, each pixel value corresponding to a pixel of 
an image representation; 
a comparing means for comparing each pixel value 
from the first image memory with a preselected 
criteria; 

a filter function altering means for selecting charac- 
teristics of a filter function for each pixel to be 
filtered in accordance with the comparison of the 
pixel value with the preselected criteria, the pixel 
value altering means being operatively connected 
with the comparing means; and, 

a first filtering means for operating on the pixel values 
of each pixel and pixels thereadjacent with the 
selectively altered filter function, the first filtering 
means being operatively connected to the first 
image memory, whereby the filter function with 
which each pixel is filtered varies in accordance 
with its pixel value. 

9. The apparatus as set forth in claim 8 further includ- 
ing: 

a second image memory for storing pixel values of a 
second image representation; 

a second filtering means for operating on the pixel 
values from the second image memory for each 
pixel and pixels thereadjacent with the selectively 
altered filter function, whereby the filter function 
with which each second image memory pixel is 
filtered varies in accordance with the comparison 
of the pixel value of the corresponding pixel of the 
first image memory with the preselected character- 
istic; 

a first filtered image memory means operatively con- 
nected with the first filtering means for receiving 
filtered first image pixel values therefrom; 

a second filtered image memory operatively con- 
nected with the second filtering means for receiv- 
ing filtered second image pixel values therefrom; 

a transform means operatively connected with the 
first and second filtered image memory means for 
transforming the filtered first and second image 
pixel values into corresponding pixel values of a 
composite image representation; and, 

a composite image memory means operatively con- 
nected with the transform means for storing the 
composite image pixel values. 

10. The apparatus as set forth in claim 9 further in- 
cluding a transform altering means for selectively alter- 
ing the transform applied by the transform means, the 
transform altering means being operatively connected 
with at least one of the first and second filtered image 
memory means for selectively altering the transform in 
accordance with the pixel values thereof. 

11. A method of diagnostic imaging comprising: 

(a) generating first and second image representations 
of the same region of a subject; 
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(b) filtering the first image representation with a first create a filtered pixel value for a corresponding 
filter function to produce a filtered first image rep- P*el of a filtered image representation; 
resentation* W rc P catm S stc P s (M through (d) for other pixels of 

(c) filtering the second image representation with the **^ g ?/ e f ^ - , « v • u 
r * <m* c ♦« „ ru^j e^^^j 5 16. The method as set forth in clam IS wherein each 
first filter taction to produce a filtered second ^ & matfix ^ wnerdn the filter 

image representation and, function selectin selectin the value$ of 

(d) transforming the filtered first and second image dement$ ^ d £ herein ^ mterin step mcludes 
representations with a first transform function into convolving pixel values of the mtigc representation 
a first composite image. 10 ^ ^ matrix. 

12. The method as set forth in claim 11 further includ- 17 ^ met h 0< i 33 set forth in claim 15 further includ- 
ing filtering the first and second image representations re peating steps (a) through (e) for a second image 
with a second filter function to produce second filter representation such that first and second filtered image 
function filtered first and second image representations representations are generated and transforming the Al- 
and transforming the second filter function filtered first tS tered image representations into a plurality of compos- 
and second image representations with a second trans- ite images. 

form function into a second composite image. 18. The method as set forth in claim 17 wherein the 

13. The method as set forth in claim 12 wherein the steps of generating first and second image representa- 
first and second image representations represent high tions include irradiating a common region of interest of 
and low energy specific radiographic images, respec- 20 a patient with radiation of first and second energy levels 
tively, and wherein the first and second transform rune- and measuring the intensity of radiation of each of the 
tions transform the filtered image representations into levels traversing the re^on of interest 
images representative of a first material and a second *f herem *f composite images include a soft tissue 

0 . . r , , , _ - basis imaze and a bone basis image, 

materud such that the first composite unage * a first W . A ^thod of medical dia^stic imaging compris- 

material specific image and the second composite image ^ °^ 

is a second material specific image. A getting fi^t second image representations 

14. The method as set forth in claim 11 wherein: of a common region of a patient; 

the first and second image representations are each ^ operating on both the first and second image rep- 
defined by a plurality of pixels, each pixel of the 30 resentations with (1) a first transform function and 
first image representation having a corresponding (2) a first filter function to create a first composite 
pixel in the second image representation which image; and, 

represents the same subregion of an imaged region (c) operating on both the first and second image rep- 

of interest, each image representation having a resentations with (1) a second transform function 

plurality of pixel values, each pixel value corre- 35 and (2) a second filter function which is different 

sponding to one of the pixels, and from the first filter function to create a second 

further including altering the first filter function in composite image, 

accordance with a characteristic of each pixel 20. The method as set forth in claim 19 wherein the 

value to be filtered such that the first filter function image representation generating step mcludes irradiat- 

is varied in accordance with the pixel value of each *> ™& the patient with radiation of a first energy level to 

Mel to be filtered generate the first unage reprcsentation and irradiating 

j — » - , c .. " ^ . . . the patient with radiation of a second energy to gener- 

15 A method of diagnostic imaging comprising: ^ . representation and wherein the 

(a) generating at least one image representation trinafem ^^tre selected such that the first com- 
which is defined by a plurality of pixels each with ^ pQsite fa mVeM spccific t0 ^ft ti ssuc md the 
an electronic pixel value; second composite image is material specific to bone. 

(b) for one pixel, comparing at least its pixel value 2 1. The method as set forth in claim 19 further includ- 
with a preselected pixel value criteria; mg adaptively altering at least one of the first filter and 

(c) in accordance with the comparison of step (b), transform functions and altering at least one of the sec- 
selecting one of a preselected plurality of filter 50 ond filter and transform functions in accordance with 
functions; characteristics of at least one of the first and second 

(d) filtering the pixel values of said one pixel and image representations. 

. adjacent pixels with the selected filter function to • ♦ * * • 

55 
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